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Pulse trains growing from modulated continuous waves (CWs) are considered, using solutions 
of the Hirota equation for solitons on a finite background. The results demonstrate that pulses 
extracted from the maximally compressed trains can propagate preserving their shape and forming 
robust arrays. The dynamics of double high-power pulse trains produced by modulated CWs in a 
model of optical fibers, including the Raman effect and other higher-order terms, is considered in 
detail too. It is demonstrated that the double trains propagate in a robust form, with frequencies 
shifted by the Raman effect. 


I. INTRODUCTION 


Recently, solitons on a finite background, generated 
by the fundamental nonlinear Schrodinger (NLS) equa¬ 
tion with self-focusing, which have been known for a 
long time, have drawn renewed interest as prototypes 
of rogue waves in the ocean [l|-[6|, optical fibers 
plasmas [T3[ and other nonlinear media. NLS solutions 
for solitons on a finite background can be identified as 
the Kuznetsov-Ma soliton (KMS) [13, [H , which is local¬ 
ized in the transverse direction, the Akhmediev breather 
(AB) 0, localized in the longitudinal direction, and the 
Peregrine solution (PS) [13 , which is a limit case of the 
KMS and AB, being localized in both the transverse and 
longitudinal directions. The spatiotemporal localization 
predicted by the PS solution, and dynamical manifes¬ 
tations of the KMS have been demonstrated in optical 
fibers [3, [T^- Splitting of the PS produced by deformed 
initial conditions, and the related modulational instabil¬ 
ity have been observed too 0[i3- The dynamics ac¬ 
counted for by the AB was first observed in Ref. [TsI 
and further analyzed in Ref. [9[ . Collisions between ABs 
0, their dynamics in optical fibers with a longitudi¬ 
nally tailored dispersion [^, and shaping of an optical 
frequency comb into a generator of rogue waves ia have 
been demonstrated too. In these experiments, the ABs 
and KMSs were excited by weakly and strongly modu¬ 
lated continuous waves (CWs), respectively. 

Generally, NLS solutions initiated by perturbations 
added to a finite background feature strong temporal 
compression of emerging peaks [23]. In optical systems, 
such dynamical scenarios can be used for the generation 
of supercontinuum, and of high-power ultrashort pulses 
and pulse trains [TsI, [234^ . 

To create high-repetition-rate and high-quality pulse 
trains, the delay-line interferometer was used to elim¬ 
inate the finite background around the AB and allow 
doubling of the repetition rate of the pulse trains [2^ . 


A possibility of spontaneous transformation of ABs into 
solitons and ultrashort pulses generated W the modula¬ 
tion instability have been considered too [^, [2^ . Extrac¬ 
tion of high-power optical pulses from the PS by means 
of the spectral-filtering method was proposed in Refs, 
[sol-ls^. However, these results are only suitable for pi¬ 
cosecond pulses with the width ^ 5 ps [13, [13 • In the 
present paper, using soliton solutions on the finite back¬ 
ground for the Hirota equation, we construct high-power 
strongly compressed pulse trains induced by modulated 
CWs in the subpicosecond regime, an appropriate model 
for which must include higher-order terms, added to the 
usual NLS equation. 

The paper is organized as follows. In Sec. H, the Hi¬ 
rota equation, which models the pulse propagation in the 
nonlinear fiber, is presented, and solutions for solitons on 
the finite background are recapitulated for this equation. 
In the same section, using AB and KMS solutions for the 
Hirota equation, high-power compressed pulse trains in¬ 
duced by modulated CWs are demonstrated in the model 
with the higher-order terms, making use of juxtaposition 
of two solutions with a half-period temporal shift and 
opposite signs. The influence of the Raman effect on the 
high-power pulse trains is considered in Sec. HI. Conclu¬ 
sions are summarized in Sec. IV. 


II. HIGH-POWER PULSE TRAINS PRODUCED 
BY THE HIROTA EQUATION 


As said above, the standard NLS equation is valid as 
the model for picosecond pulses with the width > 5ps, 
while in the subpicosecond or femtosecond regime higher- 
order terms should be included, to account for the third- 
order dispersion (TOD), self-steepening, and the Raman 
effect [^, [sj. This accordingly modified NLS equation 
is 
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Here A = A(z^T) is the slowly varying envelope of the 
electric field, T = t — zjvg = t — Piz, where t and ^ are 
the temporal variable and propagation distance, and Vg 
is the group velocity. Coefficients P 2 , 7, and Tr 

account for the second-order group-velocity dispersion 
(GVD), TOD, strength of the Kerr nonlinearity, nonlin¬ 
ear dispersion, and the Raman time-delay constant, re¬ 
spectively. Equation m does not include the fiber loss, 
assuming, as usual, that it may by compensated by gain 
[nisi. At the end of Sec. Ill, we produce results of 
the analysis of the model which explicitly includes the 
loss, the conclusion being that the difference is relatively 
small except for a decreasing of peak power. 

By means of rescaling, A{z,T) = r = 

TjTo^ and ^ with the temporal scale Tq = 

[1/^21 and dispersion length Ld = Eq. 

m is transformed into 
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where s = —^ 2 / |/^2| is +1 and —1 for the anomalous and 
normal GVD, respectively, and the other coefficients are 
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which, as well as their unsealed counterparts, are related 
to the TOD, self-steepening, nonlinear dispersion, and 
the Raman effect, respectively. In its general form, Eq. 
© is not integrable. However, if the Raman effect is 
omitted, i.e., tr = 0, and additional conditions 


0^4 = 6(^3, 0^4 + 0^5 = 0, (4) 

are satisfied, Eq. m becomes the integrable Hirota equa¬ 
tion [35I . Here we only consider the case of the anomalous 
GVD, i.e. 5 = 1. In this case, the solution of the Hirota 
equation for a soliton placed on top of a finite background 
can be compactly written as [^, 


q{C,T) 


^ ^ 2 (1 — 2a) cosh(6^) + ibsmh{b^) 
V^cos[co’(r + c^)] — cosh(6^) 
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uj = 2\/r^-^, b = c = 2(1 + 4a)(a3, (6) 

where a is a real constant parametrizing this class of the 
solutions 0- When 0^3 = 0, solution dU goes over into 
its NLS counterpart |1Q| . 

Eor 0 < a < 1/2, constants cj, b and c in Eq. ^ are 
real. In this case, the solution describes the AB, and the 
maximally compressed pulse train is attained at point 
^ = 0, in the form of 
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FIG. 1: (Color online) (a) The maximally compressed pulse 
train at z = 8.4206 km, juxtaposed with its counterpart with 
the opposite sign, delayed by a half a period. Here the blue 
dotted and green dashed-dotted curves depict, respectively, 
the real and imaginary parts of the numerical result, while 
the black solid curve shows the exact solution, which actu¬ 
ally overlaps with the blue dotted curve, (b) The intensity 
distribution of the resulting double pulse train, produced as 
the superposition of the two original trains, (c) The evolution 
of the double train with zero background, (d) The evolution 
of the peak power of the double train with zero background. 
These results correspond to a = 0.4 in Eq. (|6|) and S = 0.01 
in Eq. (|9|). 


with period 27r /uj and background power 

Pbg = (2\/^-l)2, (8) 

which corresponds to cos(ci;r) = —1. The analysis re¬ 
ported in Ref. [2^ has demonstrated that, in the range 
of 0 < a < 0.2, the evolution of the maximally com¬ 
pressed pulse train can be frozen with the help of a specif¬ 
ically tailoring the GVD profile, which leads to formation 
of a quasi-stable train of fundamental solitons, while at 
a > 0.2 a delay-line interferometer can be used to create 
a high-power pulse train with the double repetition rate 
and zero background, as shown in Ref. 

A natural question is whether such a double-repetition- 
rate high-power train can stably propagate in the frame¬ 
work of the fiber model. To address this issue, it is nec¬ 
essary to choose suitable initial conditions, which can 
excite a maximally compressed pulse array similar to the 
AB profile ([7]) . It has been shown that such a maximally 
compressed pulse train can be excited by a modulated 
cw pm, corresponding to the initial condition 

g'(0,T) = vTT^~cos~(cjr), (9) 

where 6 is a small modulation amplitude, and uo is the 
modulation frequency, taken as per Eq. 

We numerically simulated the integrable Hirota equa¬ 
tion with the initial modulated CW state (j9j), where the 
system parameters are adjusted to the SME-28 fiber with 
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= —21.4 ps^km /^s = 0.12 ps^km and 7 = 1.2 
W“^km“^ around 1550 nm ini, and the initial power 
Po = 0.7 W. The simulations are necessary, as generic 
solutions of equations cannot be obtained in an explicit 
analytical form. In fact, high-power pulse trains grow¬ 
ing from the modulated CW do not necessarily reduce 
to exact breather solutions of the Hirota equation, such 
as KMS [ 13 . In the present case, the normalized pa¬ 
rameters are Tq = 5.0474 ps, = 1.1905 km, and 
0^3 = 1.8516 X 10“^, respectively. The simulation demon¬ 
strates that the maximally compressed pulse train is at¬ 
tained at 2 ; = 8.4206 km, as shown in Fig. [D^a). At 
this position, we superpose the train and its counter¬ 
part with the opposite sign, delayed by half a period, to 
produce a double pulse train, which is displayed in Fig. 
[TJb). The procedure can be experimentally realized by 
using a delay line and a piezo-electrical device to impose 
the necessary delay of t he p ulse train and the tt phase- 
shift, as shown in Ref. [^. Then, we let the resulting 
pulse train propagate in the fiber, see the corresponding 
stable regime in Fig. [He). Unlike its exact counter¬ 
part, the numerical solution at the point of the strongest 
compression includes an imaginary part [shown by the 
green dashed-dotted curves in Fig. [Ha)], whose maxi¬ 
mum value is 0.0985 ~ 1% of the amplitude in the real 
part, which is a small perturbation that was included to 
test the stability of the propagation of the pulse train 
with the zero background. Because the soliton number, 
N = [7P(Ar/1.763)^/|/32|)]^/^, takes value 1.3617 for 
the present parameter set, which exceeds 1, individual 
pulses in the resulting train, with peak power P and full- 
width at half-maximum (FWHM) Ar = 3.3878ps, ex¬ 
hibit oscillations, as shown in Fig. ltd). We stress that 
the oscillatory regime is a dynamically stable one, as is 
clearly shown by the simulations. 

Next, we turn to the case of a > 1/2, which corre¬ 
sponds to the high background power, as per Eq. (©• 
In this case, Eq. yields imaginary uj and 6, convert¬ 
ing the hyperbolic and trigonometric functions in Eq. (|5]) 
into trigonometric and hyperbolic ones, respectively, thus 
swapping localization and periodicity features of the so¬ 
lution. It now represents a bright oscillating pulse prop¬ 
agating on top of a finite background, which is gener¬ 
ally classified as KMS. Accordingly, the maximally com¬ 
pressed pulse is attained at ^ = 0 


^max('^) — 


1 — 4a + ^/^ cosh.{ujr) 
y/^cosh{(jT) — 1 


( 10 ) 


with UJ = 2(2a — 1)^/^. This pulse has a single peak at 
T = 0 [^max(O) = — (1 + 2\/^)], and background value 
g = 1 as |r| ^ oc, which makes the solution completely 
different from the AB. 

In principle, the maximally compressed pulse can be 
created by a small localized sech-type perturbation added 
to a CW background 1^. On the other hand, it was 
demonstrated in Ref. [1^ that such a pulse can be ex¬ 



FIG. 2: (Color online) (a) The maximally compressed pulse 
train at z = 3.7080 km and its counterpart with the opposite 
sign and temporal delay tt/U. (b) The intensity distribution 
of the train produced by the juxtaposition of both original 
trains, (c) The evolution of the resulting pulse train with 
zero background, (d) The evolution of the peak power of the 
resulting pulse train with zero background. Here a = 0.6. 
The notation adopted in the panels is the same as in Fig. [T] 


cited by a strongly modulated CW in the form of 

g(0, r) = 1 + A [1 + cos (Ur)] (11) 

with U = cD cos“^[(\/4a — 2v^ + 1 — ~ 1) “ 

l]/cosh“^{l[(4a-2)/(y/4^r^272lTT- 1) - l]/V^} 

and A = \/^ — 1, which has the peak power and FWHM 
identical to those of the intensity profile, |gmin(^)P = 
{l + [2(2a — l)]/[\/^cosh(a;(r+C7r/6)) + l]}^, at the point 

of the minimal intensity, ^ = 7r/6, where b = [8a(2a — 
1)]V2, 

Starting from input m, we have simulated, similar 
to what is shown above in Fig. [H the evolution pro¬ 
duced by the juxtaposition of the maximally compressed 
pulse trains and its counterpart with the opposite sign 
and half-a-period delay. The results, summarized in Fig. 
[21 show that, for our choice of the parameters, the max¬ 
imally compressed pulse train is attained at z = 3.7080 
km [Fig. |2][a)], with individual pulses in the train be¬ 
ing almost identical to the profile of the maximally com¬ 
pressed pulse given by Eq. m- Thus, combining the 
pulse train with its counterpart with the opposite sign 
and time delay tt/U, one can construct a double pulse 
train with zero background [see Eig. Efb)]. It propa¬ 
gates stably, exhibiting oscillatory behavior of the pulses, 
as the respective soliton number is N = 1.2761, with 
EWHM Ar = 2.7927 ps, which again exceeds 1, see Eigs. 
He) and (d). Similarly, unlike the exact solution, the nu¬ 
merical result includes an imaginary part [see the green 
dash dotted curves in Eig. Efa)], with maximum value 
0.0583 ~ 0.76% of the amplitude of the real part, which 
does not affect the stable propagation of the pulse train 
with zero background. 
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Finally, it should be pointed out that, in the case of 
a 1/2, solution (|5]) of the Hirota equation reduces to 
the PS in the form of 




4 (1 + 2i0 




which is a superposition of the CW solution and a ratio¬ 
nal fraction function, forming a maximally compressed 
pulse at ^ = 0 [13, Si • This implies that the Peregrine 
rogue wave, i.e., the maximally compressed pulse, can 
be excited by a small localized (single-peak) perturba- 
tion placed on top of the CW background [^. Thus, 
the high-power pulse extracted from the Peregrine rouge 
wave can be realized by means of the spectral-filtering 
method, displaying a breather-like behavior [Siiii. 



III. THE INFLUENCE OF THE RAMAN 
EFFECT ON HIGH-POWER PULSE TRAINS 

In the previous Section, we have considered the gener¬ 
ation and propagation of the pulse trains growing from 
the CWs as solutions of the Hirota equation, i.e., Eq. 
@ under the integrability relations while the Ra¬ 
man effect was not be included {tr = 0). This model 
does not precisely apply to real optical fibers, but the 
exact solutions provided by it help one to understand 
the pulse-train dynamics in more general models. 

In this section, we report results of simulations of 
Eq. ([2]) with realistic values of parameters ([3]), includ¬ 
ing a nonzero Raman coefficient, and without impos¬ 
ing integrability conditions (|4)): 0^3 = 1.8516 x 10“^, 
0^4 = —1.6292 X 10“^, 0^5 is determined by the delayed- 
nonlinear-response parameter and tr = 5.9437x 10“^ 
for the usual value of the Raman time constant, Tr = 3 
fs. Actually, is a small quantity which has little effect 
on simulation results, therefore we set aj^ = 0. Thus, 
we focus on the influence of the Raman effect on the dy¬ 
namics of the high-power pulse trains growing from the 
modulated CWs introduced as per Eqs. m and m- 
In the analysis presented above the modulation ampli¬ 
tude A and frequencies uj and U in Eqs. m and (ED 
were taken as functions of parameter a, making it possi¬ 
ble to compare the results with the exact solutions pro¬ 
duced by the Hirota equation, as shown in Eig. [T]^a) and 
Eig. [2fa). However, in the realistic fiber-optic model, 
the modulation amplitudes and frequencies can take ar¬ 
bitrary values. Here, we use inputs given by Eqs. 
and m with different values of the parameters, aiming 
to excite the maximally compressed pulse trains, and per¬ 
form subsequent simulations of the corresponding double 
pulse trains in the framework of Eq. ©• 

Eig. [ 3 ] shows the evolution of the train initiated by 
input ([9]) with 6 = 0.03 and 00 = 1. In Eig. [3l^a) 
one can see that the modulated CW initially undergoes 
strong temporal compression and increase of the peak 
power, attaining the maximally compressed state of the 
train at z = 6.5277 km. At that stage, the double 


FIG. 3: (Color online) (a) Power profiles at different propaga¬ 
tion distances for the pulse train generated by the modulated 
CW given by Eq. (|9|). (b) The corresponding power spectra, 
(c) The evolution of the peak power. Here the parameters are 
S = 0.03, UJ = 1, Tr = 3 fs, and = 0. 



FIG. 4: (Color online) The same as in Fig. [3l but for the 
pulse train generated by the input in the form of Eq. (EJ. 
The parameters are A = 0.13, Q = I, Tr = 3 fs, and = 0. 


pulse train is constructed, using the delay-line interfer¬ 
ometer (the method is outlined in the previous Section, 
see also Ref. il), which robustly propagates along the 
fiber. Eig. [H^b) displays the concomitant evolution of 
the pulse-train’s power spectrum, showing that the spec¬ 
trum initially expands until the point of the maximum 
compression, at which the spectral distribution is reduced 
to a half of the original spectrum • Subsequently, the 
spectral width remains essentially constant, while the fre¬ 
quencies feature the red shift induced, as usual, by the 
Raman effect [ 31 , 133 . Eor the presently chosen parame¬ 
ters, the frequency shift of the train is about 1.1967 THz 
at 60 km. In addition, the evolution of the peak power. 
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FIG. 5: (Color online) (a) The position of the maximally 
compressed pulse train, (b) the peak power, (c) the FWHM 
of the central pulse in the double pulse train, and (d) the 
Raman-induced frequency shift after passing 60 km in the 
fiber, versus the modulation intensity 5 or A^ with the black 
solid and blue dotted curves pertaining to inputs ([9]) and CB, 
respectively. Here the parameters are Q = cj = 1, Tr = 3 fs, 
and ai, — 0. 


which exhibits an oscillatory behavior, is displayed in Fig. 

He). 

Fig. m displays the evolution initiated by the input 
given by Eq. m with = 0.13 and ft = 1. This sce¬ 
nario exhibits three differences in comparison with Fig. 
[3l One is that the maximally compressed pulse train 
is attained at a shorter distance, i.e., 2 ) = 2.9655 km, as 
seen in Fig. lU^a). The other is that the initial spectrum is 
narrower, while the Raman-induced red frequency shift is 
larger, as seen in Fig. lU^b), in which the frequency shift 
of the train is about 2.3854 THz at 60 km. The third 
difference is a higher peak power, as seen in Fig. SJ^c). 

To summarize the dynamical results outlined above, 
in Fig. [5] we present the largest-compression position Zc 
for the pulse trains, their peak power, the FWHM of the 
central pulse in the double pulse trains constructed with 
the help of the delay-line interferometer method, and the 
corresponding frequency shift accumulated after the 60 
km long propagation, as functions of parameters S and 
A in inputs (j9j) and m, respectively. 

In Fig. [5fa), one can see that, for the same modulation 
intensity, the position of the strongest compression corre¬ 
sponding to initial condition m is always smaller than 
its counterpart corresponding to input ^ . Interestingly, 
for the double pulse train induced by the latter input, 
the peak power and the FWHM of the central pulse are 
almost constant for the present parameters, which leads 
to the nearly constant frequency shift at 60 km. On 
the other hand, for the double pulse train induced by 
input m, the peak power increases, and the FWHM 
decreases, with the increase of the modulation intensity, 
which implies that this input can be used to generate 
the pulse train with the higher peak power and narrower 


FIG. 6: (Golor online) The same as in Fig. [4| except for 
A = 0.55. 

width. However, in the same case, the frequency shift 
at 60 km increases with the increase of the modulation 
intensity in a staircase fashion, as shown in Fig. [5fd). 
Furthermore, from Figs. Hb) and (c) one can find that 
the corresponding soliton numbers are in the range of 
1.2592 < N < 1.3412, each pulse in the double pulse 
trains eventually reshaping into the fundamental soliton, 
as seen in Figs. ^c) and He). Comparing with Figs. 
[Tfd) and[2fd), one can see that Raman effect is helpful 
for the formation of solitons over a shorter propagation 
distance, as was shown in a different context in Ref. [39|. 

From Fig. [5fc) it can be seen that the pulse width 
in the trains is less than 2 ps (here the pulse width is 
taken equal to FWHM/1.763 [s^), and it may fall in 
the subpicosecond range for larger values of modulation 
intensity A. These widths are small enough to make it 
necessary using the extended model based on Eq. o, 
rather than the simple NFS equation which does not in¬ 
clude higher-order terms. To additionally illustrate this 
point, we have performed simulations of the propaga¬ 
tion dynamics of the train initiated by input m with 
A = 0.55 in the subpicosecond range. In this case, the 
maximally compressed pulse train is attained at a shorter 
distance ^ = 1.0114 km in comparison with Eig. 01 Sub¬ 
sequently, the use of the delay-line interferometer pro¬ 
duces the double pulse train with the EWHM of the cen¬ 
tral pulse 1.4083 ps (the corresponding pulse width is 799 
fs) and peak power 26.3165 W. The corresponding robust 
evolution of the double pulse train is shown in Eig. [Sfa). 
Comparing with the results shown in Eig. 01 we conclude 
that the spectrum range is wider, the frequency shift is 
larger, and the peak power is higher because the pulse 
width becomes narrower, as seen in Eigs. [D^b) and (c). 

Note that the above considerations disregarded fiber 
loss, assuming that it might be compensated by gain. 
Here, we display an example which explicitly takes the 
loss into account, and the results are shown in Eig. [71 in 
which the maximally compressed pulse train is attained 
at a somewhat longer distance, i.e., z = 3.1306 km (it 
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FIG. 7: (Color online) The same as in Fig. IH but with fiber 
loss 0.19 dB/km. 


(Akhmediev breather) and KMS (Kuznetsov-Ma solu¬ 
tion), the double high-power pulse trains with zero- 
background have been constructed. The results demon¬ 
strate that the double trains, built as superpositions of 
the maximally compressed pulse trains with opposite 
signs and a half-period time delay, propagate preserv¬ 
ing their shape. The dynamics of the trains generated by 
the modulated CWs in the realistic fiber model has been 
simulated too, demonstrating that they can propagate ro¬ 
bustly, with the red frequency shift caused by the Raman 
effect. Comparing different initial conditions, we have 
concluded that input m is more suitable for generating 
the pulse trains with a higher peak power and narrower 
width, the corresponding Raman-induced frequency shift 
being larger too, exhibiting a staircase dependence on 
the CW-modulation parameters. The settings analyzed 
in this work suggest an experimental realization in non¬ 
linear optical fibers. 


was 2.9655km without fiber loss) and the peak power 
of the corresponding double pulse train is decreased to 
9.6711 W, from 11.6262 W in the lossless fiber. Com¬ 
paring with Fig. m we see that the propagation scenario 
is not strongly affected by the loss. Of course, the loss 
produces stronger effects over large propagation distances 
(> 30 km), in which case periodic compensation provided 
by gain should be explicitly taken into account. 
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